The Berry phase of dislocations in graphene and valley conserving decoherence 
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We demonstrate that dislocations in the graphene lattice give rise to electron Berry phases equiv- 
alent to quantized values {0, ±|} in units of the flux quantum, but with an opposite sign for the 
two valleys. An elementary scale consideration of a graphene Aharonov-Bohm ring equipped with 
valley fllters on both terminals, encircling a dislocation, says that in the regime where the intervalley 
mean free path is large compared to the intravalley phase coherence length, such that the valley 
quantum numbers can be regarded as conserved on the relevant scale, the coherent valley-polarized 
currents sensitive to the topological phases have to traverse the device many times before both 
valleys contribute, and this is not possible at intermediate temperatures where the latter length 
becomes of order of the device size, thus leading to an apparent violation of the basic law of linear 
transport that magnetoconductance is even in the applied flux. We discuss this discrepancy in the 
Feynman path picture of dephasing, when addressing the transition from quantum to classical dissi- 
pative transport. We also investigate this device in the scattering matrix formalism, accounting for 
the effects of decoherence by the Biittiker dephasing voltage probe type model which conserves the 
valleys, where the magnetoconductance remains even in the flux, also when different decoherence 
times are allowed for the individual, time reversal connected, valleys. 

PACS numbers: 72.10.Fk,73.23.-b,73.63.-b,03.65.Yz 



I. INTRODUCTION 

Since electrical conductance is even under time rever- 
sal, it has to be that magnetoconductance is an even func- 
tion of the applied magnetic field that breaks time rever- 
sal invariance. This elementary Casimir-Onsager rela- 
tion requires equilibrium conditions such that the trans- 
port is in the linear response regimeJ^i^i^ 

Here we will present an example suggesting that in 
the case of finite temperature quantum transport, linear 
response might run into a singular limit: although the ex- 
ternal conditions are perfectly within linear response, the 
parts of the current that are governed by quantum me- 
chanics cannot equilibrate in a true sense because some 
quantum number is effectively conserved, with the net ef- 
fect that these coherent currents feel an "arrow of time" 
negating the Onsager relations associated with true equi- 
librium. This might be a more general truth, but we will 
limit ourselves here to the specific case of graphene where 
we have to employ a whole array of properties specific to 
graphene, to come up with a design that might exhibit 
the aforementioned effect. We stumbled on this story in 
trying to find out how to turn a topological phase, that 
is most natural to graphene Dirac electrons, into an ob- 
servable quantity. 

The effect of topology on electronic properties is tied to 
the topological features of the underlying atomic lattice. 
These are the dislocations and disclinations. Although 
disclinations are rather unnatural according to the stan- 
dard theory of elasticity (or plasticity),^ the global in- 
fluences they exert on graphene's Dirac electrons^'^ have 
been relatively thoroughly studied, with a special focus 
on the similarities with the holonomy structure of funda- 
mental Dirac electrons in a curved space-t ime .— 10,11,12 

However, dislocations have been largely ignored ; ^^1^^ 
although these are ubiquitous topological defects in any 



solid. In contrast to disclinations they require only finite 
energies to be created, so that it is virtually impossi- 
ble to prepare a crystal that contains no dislocations at 
all. These have not been found in the graphene flakes 
produced by the Manchester method, likely because 
dislocated graphene does not survive this rather violent 
method of preparing a sample. With more sophisticated 
manufacturing methods it is expected that graphene dis- 
locations will be abundant 

A dislocation, due to its topological nature, exerts in- 
fluence also far away from the core. The question arises 
as to what happens to a quantum-coherent graphene 
Dirac electron that is transported around it. We ana- 
lyze this problem in Section [TTl the outcome is a holon- 
omy structure of pleasing simplicity. The dislocation is 
the topological defect associated with translations^^ and 
since translations are Abelian, the holonomy is akin to 
the holonomy associated with electromagnetism — the 
Aharonov — Bohm (AB) phase. A crucial difference is 
that the Dirac electrons feel a (pseudo-)flux of the same 
magnitude, but opposite sign in the two valleys, which 
is a consequence of dislocations leaving the system time 
reversal invariant. In addition, the topological charge of 
the dislocation appears in quantized units corresponding 
to fractions {0, of the magnetic flux quantum. 

For detection purposes, one envisages a typical AB ex- 
periment where the dislocation is put in the middle of a 
ring (Fig. [1]). The AB oscillations are infiuenced by the 
presence of the dislocation holonomy, and we will dis- 
cuss this in Section HVl It turns out that the dislocation 
topological phase could in principle be measured, after it 
is disentangled from the elastic scattering of impurities 
by disorder averaging. Its effect is connected to the AB 
oscillation amplitudes, which are in practice less reliable 
due to the standard mesoscopic clutter of the oscillations. 

If the current was carried exclusively by electrons in 
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one vallejsi^ the situation would be quite different, since 
these sense the dislocation Berry phase as indistinguish- 
able from a real magnetic flux. Abstractly, it seems the 
dislocation Berry phase could thus cause the offset of the 
magneto-oscillations, which would violate the Onsager 
relation. We therefore consider the concrete possibility 
of valley filters installed at the input and output termi- 
nals of our dislocated AB ring (Fig.[T]). The time reversal 
invariance puts a constraint on the general workings of 
the valley filter: when it is perfectly transparent for elec- 
trons in valley coming from the left (thus completely 
reflecting the K- valley) , a electron impinging on it 
from the right will be unitarily backscattered to the K 
valley. Deep in the quantum regime where the phase co- 
herence length is large compared to the size of the ring L, 
long Feynman paths encircle the ring many times, hav- 
ing ample opportunity to explore the "backside" of the 
valley filters, with the effect that the quantum current 
equilibrates over the two valleys, restoring the evenness 
of the magnetoconductance. When temperature rises, 
the phase coherence length shrinks and becomes of order 
of L. The coherent part of the current that is sensitive to 
the topological phases can still be detected but now it is 
dominated by Feynman paths that traverse the ring only 
once. These can no longer explore the backside of the 
valley polarizers, and so can no longer sense that time 
reversal is unbroken, with the consequence that the mag- 
netoconductance becomes uneven. We will address this 
more quantitatively in Section IV Al The simple essence 
of the argument is the observation that even in a lin- 
ear response measurement, the quantum coherent part 
of the current cannot reach a true equilibrium. The un- 
derlying assumption is that the electrical currents are 
conserved separately for the two valleys everywhere in 
the device, except at the valley polarizers. Since these 
are separated in space by a length L, this current can be 
regarded as effectively conserved when the phase coher- 
ence length becomes of the order of L for the purposes of 
quantum coherent phenomena that depend on the con- 
servation of valley current. This quantum conserved cur- 
rent acts in analogy with the role of conservation laws in 
conventional hydrodynamics to prohibit the system from 
reaching equilibrium. 

The argument as presented implicitly rests on the lan- 
guage of Feynman paths and there are precedents known 
where qualitative arguments of that kind can be quite 
misleading with regard to quantum transport A su- 
perb theory describing transport deep in the quantum 
regime is the Landauer-Biittiker scattering matrix for- 
malism and we will address the workings of our device 
in this language in Section IV Bl It seems that the for- 
malism is inherently static, revolving around elastic scat- 
tering which is sufRcient at zero temperature, but at 
finite temperature the role of imaginary time becomes 
central in properly accounting for the effects of inelastic 
scattering. Among the various attempts fi^i^S the volt- 
age probe approach to incorporating dephasing2ii22 sug- 
gested by Blittiker is particularly prominent. It amounts 




FIG. 1: The modified graphene Aharonov-Bohm device. This 
is the usual ring pierced by an external magnetic flux $, but 
now with a dislocation with Burgers vector b in the center 
acting as a pseudoflux on the electrons with a definite valley 
number. Both leads are equipped with a valley polarizer: ide- 
ally these transmit fully, say, electrons in the K+ mode (solid 
lines) moving from left to right, while K- electrons (dashed 
lines) propagating in the same direction are refiected to K+ 
mode moving in the opposite direction, as required by time 
reversal invariance. In Sec. IV Al it is argued that at tempera- 
tures such that the device size is of order of the phase coher- 
ence length, only Feynman paths traversing the device once 
contribute to the magneto-oscillations: K+ modes moving 
from left to right (solid lines) sense the direction of the Burg- 
ers vector in a way that is opposite to the K- modes moving 
from right to left (dashed lines), and this implies that the dis- 
location pseudofiux offsets the magneto-oscillations. At low 
temperatures the long Feynman paths explore the "backside" 
of the polarizers and Biittiker's law is restored. 



to attaching an extra terminal to the coherent quantum 
device, with the effect of scrambling the phase of the 
waves entering this phantom reservoir. This has a re- 
spectable track record with regard to correctly modeling 
the effects of dccoherence on quantum transport (e.g., 
Ref. 18; 22, 23, 24, 25). We straightforwardly extend this 
method to the present device by requiring that the de- 
phasing reservoirs do not affect the valley quantum num- 
ber, assuming the intervalley inelastic time to be in- 
finitely long. As long as time reversal and unitarity of 
scattering are present, it follows generally from this for- 
malism that magnetoconductance is even^ a fact in this 
context referred to as Biittiker's theorem. We prove that 
this holds even when different dephasing times are al- 
lowed for the two valleys, which are connected by time 
reversal. Furthermore, we explicate how the dislocation 
phase signature in the AB oscillations remains the same 
as in the zero temperature calculation. 

We hope that this story will motivate experimentalists 
to realize our device in the laboratory. It appears to us 
that the matters at stake cannot be decided by theoret- 
ical means alone, as we will substantiate in the rest of 
this paper. 
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II. ELECTRON BERRY PHASE AND THE 
BURGERS VECTOR OF DISLOCATIONS 

One of the two possible topological crystal defects, the 
dislocations, are omnipresent in crystals in general. A 
dislocation is in principle obtained by the Volterra con- 
struction as follows: a semi-infinite strip of unit cells 
is removed from a crystal, and the open edges are glued 
back together along the Volterra cut, leaving some imper- 
fections at the original beginning of the strip (the core), 
see Fig. [2] Tracing a closed loop around the defect core, 
but drawing it in the perfect lattice, one finds a non- 
closure, equal to some lattice vector — the Burgers vec- 
tor. This persists for loops of arbitrary size, and so the 
effect of the defect on electron wavefunctions is global and 
long-ranged. This property enables one to model the de- 
fect as a nontrivial boundary condition on the wavefunc- 
tion at the Volterra cut, which can be imposed by a gauge 
field in a reversal of the usual argumentation for appear- 
ance of the Aharonov-Bohm effecti^ The difference with 
the case of disclinations ,—1^121 the other topological crys- 
tal defect, in which one cuts out a pie segment of the 
lattice, is that instead of rotating the electron spinor, 
under the influence of the translational dislocation,^ the 
spinor is translated by the Burgers vector to maintain 
single- valuedness. 

It is known that disclinations cause a deficit angle in 
loops circling the core, which in graphene can be any of 
the five multiples of ±-|-, producing a variety of physical 
effects)^ and are interesting primarily because of their 
occurrence in nanocones and fuUerenes.— 

The theoretical study of dislocations, however, has so 
far been scarce. Random distributions of dislocations 
have been discussed from the perspective of their statisti- 
cal influence on coherence and electron propagation! ^'^i 
We will here address a different set of phenomena associ- 
ated with their topology. We will show that although the 
topological charge of a dislocation could be any lattice 
vector, they act as a simple Aharonov-Bohm flux located 
at the defect core, of opposite signs in two valleys. They 
fall into three possible classes — a trivial one (zero fiux), 
and two of opposite sign (±i flux). 

Let us start by reviewing the standard low-energy, con- 
tinuum description of the graphene electron states com- 
ing from the Pz carbon orbitals i^^i^^'^^i'^^ The two "val- 
ley" Dirac points are labeled by K± = ±K (Fig. (2]), and 
the unit cell contains two atoms (labeled A and B) , yield- 
ing a total of four massless states. In this basis the wave- 
functions are described by a slowly-varying four compo- 
nent spinor. Operators acting on the A and B states 
without mixing the K± valleys are written as Pauli ma- 
trices (Ta, a G {1:2,3}; while the valley degeneracy is 
tracked by a second set of r Pauli matrices. To lowest 
order this yields the usual Dirac Hamiltonian, 

H = -i[{K-d)T3(g)ai + {A-d)l(g}a2], (1) 
where the energy is measured in units of hvp, K_ is the 



normalized K vector and A the normalized vector con- 
necting the A and B sites (see Fig. [2]). 




FIG. 2: The electronic structure and dislocations in graphene. 
By removing rows of unit cells a dislocation with Burgers vec- 
tor b is created. The ellipses indicate which unit cells can be 
removed to obtain a "trivial" dislocation not carrying a net 
topological charge, as can be seen for instance by counting 
the phases of the Bloch waves. An arbitrary Burgers vector 
starts from the central square and reaches the center of some 
hexagon, and this labels the dislocation's class: bold hexagon 
sides represent trivial dislocations (d = 0), grey shade the 
d = I class, and white fill the d = — i class. Graphene's 
Dirac electrons carry unit cell [A/B) and "valley" K± in- 
dices. The phases of the Bloch waves of the states on 
the rows of the defect-free lattice are indicated at the right in 
terms of z = ea;p[j27r/3]. By creating the Volterra cut associ- 
ated with the dislocations it follows that the Dirac electrons 
experience topological phase jumps of ^ and — ^ for dislo- 
cation class 1/3 and —1/3, respectively. The K- states expe- 
rience the opposite phase jump. Note that the phase jumps 
are independent of the A/B quantum numbers, because the 
dislocation does not affect the intra unit-cell structure. 

We now consider the influence of dislocations on such 
Dirac fermions, associated with the translation by a 
Burgers vector b at the modified boundary arising from 
the Volterra cut. The components of are coefficients 
multiplying the Fermi states, K±A/B, being Bloch 
eigenstates of the crystal lattice, and a translation by 
a lattice vector is therefore equivalent to a multiplication 
by the corresponding phase factor exp{iK± ■ b). This 
yields the U{1) holonomy 

U{b) ^ e'^^-''^'"' e*ir(''i-''2)r3^ (2) 

where bi and 62 are the integer components of the Burg- 
ers vector b in the lattice basis (see Fig. [2]). The dislo- 
cations thus separate into three equivalence classes, la- 
beled by d e {0, i,— |}, with 3d = (61 — 62) mod 3, 
where the period of 3 follows from the periodicity of the 
Fermi states (see Fig. [2]) . Different from the case of discli- 
nations,^'^ this is independent from the A/B sublatitce 
pseudospin quantum number since translations carry no 
information on the structure inside the unit cell. Instead, 
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this phase does depend on the valley quantum number 
in a simple way: the absolute magnitude is the same and 
the phases in the two valleys just differ by a minus sign. 

Avoiding the dislocation core (which shrinks to a point 
in the continuum limit), its influence can be encoded by 
adding a U{1) gauge coupling to the Dirac Hamiltonian 
in Eq. (P), 

Hd^si^H-i^-^{K-b)T3, (3) 
Zirr 

where r and ip are the standard polar coordinates, tak- 
ing the dislocation core as origin. The induced gauge 
field is in in precise correspondence with the one of 
an Aharonov-Bohm solenoid with flux in units of 
e/h, for the ±K valley electrons. Numerical simula- 
tions have already hinted that dislocations behave as 
pseudo-magnetic fluxes, in that they create vortex cur- 
rents around their corei^^ 

We close this subsection by discussing the role of time 
reversal invariance, being an important issue in this pa- 
per. Real magnetic fields break time reversal, expressed 
through an antiunitary operator (T) which involves com- 
plex conjugation (operator C). Time reversal applied to 
the graphene Dirac electrons exchanges the Fermi points 
and the corresponding modes in the leads. The time re- 
versal operator can be chosen as simply T = tiC. One 
has [T, U{b)] — 0, as well as [T, H] = 0: time reversal 
amounts to flipping the external magnetic field, reversing 
the direction of motion of electrons, and switching them 
to the opposite valley, while keeping b unchanged. After 
all, the lattice defect is just a complicated rearrangement 
in the lattice potential, and cannot break time reversal 
symmetry. The time reversal symmetry also dictates that 
the dislocation pseudoflux has to be of opposite sign for 
the two valleys. 



III. GENERAL PROPERTIES OF A 
DISLOCATED AB RING 

In the previous section we have shown that disloca- 
tions correspond to quantized magnetic fluxes, carrying 
however opposite signs with respect to the two valleys. 
The standard way to measure such fluxes is by measur- 
ing the conductance of an Aharonov-Bohm ring device, 
as indicated in Fig. [S) Besides the usual magnetic flux 
that can be pierced through the ring, we consider one or 
more dislocations located inside the ring. The electrons 
do not explore the dislocation cores and only communi- 
cate with the "lines of missing atoms" attached to the 
dislocation cores that cross the ring "somewhere" (the 
choice of this missing line is actually a gauge freedom on 
its own^). 

With regard to the feasibility of realizing this device 
in the laboratory, we already argued in the introduction 
that dislocations should be plentiful in graphene that is 
produced with non violent methods. Graphene struc- 
tures with a size < Ifim have been manufactured, and 



show quantum transport phenomena, including AB os- 
cillations. "^S, Concerning the final important ingredient, 
the valley polarizers, it has been suggested that the val- 
ley polarized currents could be generated using valley fil- 
ters constructed from thin strips of graphene with zigzag 
edges 

At this point one may ask how realistic is it to as- 
sume that valley currents are conserved at the mentioned 
length scales. The first issue is that the intravalley in- 
elastic scattering time should be, at a given temperature, 
much smaller than the intervalley inelastic scattering 
time, to satisfy the requirement that the intravalley phase 
coherence length becomes quite small while the valley po- 
larization is not destroyed at this temperature. The ori- 
gin of these inelastic times is of course not mysterious: it 
is rooted in Fermi-liquid electron-electron and electron- 
phonon scattering. Although we are not aware of un- 
ambiguous experimental informatiou j'^^'^^i'^^'^^i'^^'^^ it is 
widely believed that the intervalley inelastic time is in- 
deed much longer, because of the kinematical bottleneck 
that is active both for electron-electron and electron- 
phonon scattering in the form of the large momentum 
that has to be absorbed when the on-shell electrons are 
scattered between valleys. In fact, the elastic interval- 
ley scattering is more worrisome since valley quantum 
number is quite fragile, being rooted eventually in lattice 
potentials, and one expects it to be very sensitive to the 
imperfections of real life devices. 

There are indications from theoretical studies that the 
problems are manageable as long as one does not make 
the structures too narrow. The boundaries do not seem 
to play a critical role,**''*^ and there is numerical evidence 
for valley conservation in the ring geometry. Eventu- 
ally, one can contemplate even smooth terminations us- 
ing mass confinement due to potentials;^ which auto- 
matically preserve the valley. 



IV. DISLOCATED AHARONOV-BOHM RING 
AT ZERO TEMPERATURE 

The focus of this section will be on the fully coher- 
ent quantum transport at zero temperature and in this 
regime the valley filters do not have a decisive influence 
on the conductance. As announced in the introduction 
this might be different at flnite temperatures. The con- 
clusion of this section will be that when the valley cur- 
rents are conserved the dislocation Berry phase is observ- 
able in principle, but harder in practice: after inserting 
a dislocation in the ring, keeping it the same otherwise, 
especially with regard to point disorder, its presence can 
be deduced in principle from changes in the amplitude 
of magnetoconductance oscillations. When the interval- 
ley scattering length becomes smaller than L (ring arm 
length) , the electron transport carries no information any 
longer pertaining to the presence or absence of the dislo- 
cation(s). 

Let us focus on an ideal device which has ballistic 



5 



transport in the arms, the magnetic field, and the dis- 
location. The total topological phase contribution to 
the wavefunction on traversing the ring is just the sum 
of the electromagnetic ($, in units of h/e) and defect, 
Eq.([2]), pseudofluxes, since both electromagnetism and 
dislocations are governed by Abelian symmetries (C/(l) 
and translations, respectively). Starting with the case 
when intervalley scattering is assumed to be absent, while 
the valley filters of Fig. [1] are switched off, the current is 
due equally to carriers from both valleys. We learn from 
Eq. that for the nontrivial dislocations the magne- 
toconductance curve G(<I>) is shifted by ^ for carriers 
at one Dirac point, and by — ^ at the other, while the 
signs reverse on switching the dislocation class. Adding 
the two currents, each with the associated phase shift, re- 
sults in the magneto conductance G{<^+^-) + G{^~'^). 
Fourier expanding this as G($) = G(°) -l-G'^^) cos($) + . . . 
shows that the harmonics of order 3n, for n G Z, do not 
change, and all others are multiplied by a factor — i . In 
particular, the fundamental frequency oscillation (with 
period -) is halved in amplitude. This means that the 
influence of the dislocation Berry phase is quantitative, 
affecting only the amplitudes of the Fourier components 
of the AB oscillations. But these are also affected by 
point disorder, which gives rise to the standard sample- 
to-sample mesoscopic fluctuation. 

Let us address these matters quantitatively using 
the Landauer-Buttiker scattering matrix formalismi^ We 
employ a model where the polarizers and two arms of the 
ring are described by a single scatterer each, co mpletely 
analogously to the normal metal ring case in Ref. [4j!49| . 
The modes are labeled by transversal momentum and val- 
ley, while the electrons can propagate in both directions, 
inside both the left and right lead4^ The amplitudes 
of outgoing modes, O = ( ) ' ^'^'^ incoming modes, 
I = ( ) , are connected by a scattering matrix S_, with 
O = SI. The important submatrices of S_ are t and 
t., given by Or = IIl and Ol = LIr, where Il/Ir 
are columns of amplitudes of incoming (into a scatterer) 
modes from the left /right, and Ol/Or are columns of 
amplitudes of outgoing modes from the left /right, see 
Fig. IHb); thus t and ti are M x M matrices, where 
M is the number of modes in one lead. We employ 
the usual simplification of using only a single transver- 
sal mode (M = 2) for simplicity, with II = ( j- ) ' 

etc., with the expectation that the salient features of this 
model survive in the realistic case of graphene with more 
modes. In the remainder the ±K modes are labelled 
by a G {-|-, — }, and we follow the convention that the 
scattering matrices are defined by organizing the ampli- 
tudes in columns as described above. Notice that since 
K = time reversal exchanges the two valleys, 

connecting e.g. incoming (left moving) electron ampli- 
tude in one valley, to the outgoing (right moving) am- 
plitude in the opposite valley, on the same side of the 
scatterer. 

The scattering matrices used to calculate the total S_, 



are as follows: the splitter (circle in Fig. [Blja)) has a 
perfect transmission and divides the amplitude equally 
between the two ring arms, corresponding to the leads 
strongly coupled to the ring, i.e. e = in Ref. '45'; the 
scattering in ring arms (squares in Fig.[6Ka)) provides the 
necessary total flux phase upon encircling the ring, i.e. 
traversing both arms. We present the ballistic case for 
the upper arm: 

/ ae"* \ 

_ ae"* te"^*-'') 

S^ -e ^^-^ni^+d) ag-*^* 

(4) 

with t = a/1 — 7^, a = ij with 7 G [0,1], and is 
an effective phase encoding for the point disorder. The 
probability of transmission in the same valley, \t\'^, and 
the probability of transmission with scattering to the op- 
posite valley, |ap = 7^ are parametrized by 7, whose 
value corresponds to infinite intervalley scattering time 
for propagating through the arms. For the lower arm 
we then take — S_^{^ —^,d — d), as trav- 
eling from left to right must give opposite phase con- 
tributions in the two arms. The magnetoconductance 
curve is then calculated by the Landauer-Biittiker for- 
mula'* G($) = Tr{t{^)t}{^)), where t belongs to the to- 
tal scattering matrix of the device, obtained by combin- 
ing the ingredients we listed above. The matrix elements 
of the i($) matrix that determine the magnetoconduc- 
tance are obtained by explicitly solving for the outgoing 
amplitudes in the right terminal of the device, after fixing 
the incoming amplitudes in the left terminal to 1. Let 
us finally explicate some symmetry constraints on the 
scattering matrices. The unitarity of scattering implies 
expressing that particle current is conserved. 
Time reversal plays an important role in what follows, 
and it implies that for any matrix S_ (for a certain choice 
of phase relation between incoming and outgoing modes) 

/O 1 0\ 

m)^XS^{-^)X, ^= J J , (5) 

\0 1 0/ 

where the matrix X exchanges the valleys. Valleys act 
in a similar way as spins, with the spin up and down 
modes behaving similarly under time reversal. We will 
come back to this issue in the concluding section. 

Let us now discuss the characteristic features of the 
experimentally observable conductance G(<I>). The inter- 
valley scattering 7 is the important parameter and we 
first analyze the case when it vanishes. This corresponds 
to the case of Eq. ^ after setting 7 = 0. It is obvious 
that the dislocation pscudoflux just adds to the magnetic 
flux. Furthermore, the two valleys are decoupled in the 
whole device, implying that the two currents can just be 
added. We can then repeat the simple argument from 
the beginning of this section to obtain the "halving of 
amplitudes" rule. This is independent of the particular 
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point scatterer distribution, parametrized by the phase 4> 
in (HI). In Figure [3{a) we show the magnetoconductance 
without- (thick sohd hne) and with (red dashed hne) a 
non-trivial dislocation present in the ring, with one fixed 
disorder phase (p ~ 2.3, where one immediately discerns 
the main effect of the dislocation: the fundamental har- 
monic is multiplied by a factor — i. 

This example however hides a problem. Namely, a ring 
in the absence of dislocations, with a fixed disorder real- 
ization (0 = 01 ), and a dislocated ring with a different 
disorder realization (0 = (black solid and blue dot- 
dashed lines, respectively, in Fig.[3Ub)), produce different 
outcomes, and it becomes impossible to recognize a re- 
lationship between the two. The problem is that point 
disorder by itself can change the harmonic content of the 
AB oscillations in arbitrary ways. This has the effect 
that the specific information associated with the pres- 
ence of the dislocation becomes completely hidden for 
the experimentalist, who has to produce a new sample 
to compare a dislocated- with a non-dislocated AB ring, 
thereby changing the disorder configuration. 

However, the simple rule of halving the amplitude, de- 
scribed above, is rooted in topology, and it does survive 
when the point disorder is averaged over, which is a pro- 
cedure that can be implemented in practice. This fact is 
demonstrated in Fig. |4l where we show the results for the 
amplitudes of conductance harmonics, obtained after an 
averaging over the disorder phase 0; the first and second 
harmonic amplitude of the dislocated ring (red star) have 
half the value (and opposite sign) compared to the ones 
of the ideal ring (black square). 

The effect of intervalley scattering can be studied by 
switching on the a parameter in Eq. As an illustra- 
tion we show by thin gray lines in Fig. [2Ka) the change 
in magnetoconductance as we gradually decrease the in- 
tervalley scattering length; it interpolates between the 
outcomes of the ring with- and without- the disloca- 
tion. In Fig. [5] we show the evolution of the disorder 
phase averaged Fourier components, and these examples 
make it immediately clear that as the intervalley scatter- 
ing length becomes smaller than the ring size, informa- 
tion regarding the presence of the dislocation is wiped 
out completely. The physical reason is simple. Consider 
again the Feynman paths; the quantum conductance is 
governed by paths that encircle the ring many times, and 
such a long path will cross the dislocation "Dirac string" 
many times. But when the intervalley scattering length 
is short it will randomly carry a or valley iden- 
tity when it crosses the Dirac string, thereby picking up 
randomly the plus and minus dislocation Berry phase, 
with the obvious outcome that the net phase will aver- 
age away, and this means in turn that the current will lose 
all information regarding the presence of the dislocation. 

Finally, what is the specific effect of adding valley fil- 
ters at the leads in the quantum regime? The scatter- 
ing matrix describing the filter (half black rectangle in 
Fig. El^a)) perfectly transmits a — + modes from left to 




FIG. 3; Intervalley and disorder scattering dependent mag- 
netoconductance at zero temperature, (a) Consider a fixed 
disorder configuration (0 — 2.3): at infinite intervalley scat- 
tering time (a = in Eq. @), we obtain the thick black line 
in the absence of a dislocation, and the thick dashed red line 
in presence of a d = 1/3 dislocation. Notice the amplitude 
relation described in Sec. IIVI The thin gray lines show the 
evolution of the dislocated case with shortening intervalley 
scattering time. In the limit of maximal scattering (a — > 1), 
the thick black line is reached, as if no dislocation is present, 
(b) Illustration of the influence of point disorder, with no 
intervalley scattering. The solid black line (absence of dislo- 
cation) and the dashed red line (in presence of dislocation), 
at a fixed disorder configuration — 2.3, are identical to the 
ones in part (a). The dot-dashed blue line is obtained in the 
presence of a dislocation, but with the disorder configuration 
changed to = 0.1. In contrast to the case of the red line, 
the blue line has no "halving the first harmonic amplitude" 
relationship (see text) to the black line, as different disorder 
configurations can produce dramatically different AB oscilla- 
tions. 



right, and so time reversal symmetry fixes the form 



^pol 



/O 1 0\ 

1 

10 

Vo 1 o/ 



We already emphasized in the introduction that time re- 
versal symmetry implies that the backside of a perfect 
valley polarizer acts like a perfect intervalley scatterer, 
as is further illustrated in the inset of Fig. [TJ Even in 
the absence of any other source of intervalley scattering, 
this implies that valley currents are no longer conserved, 
since the long Feynman paths will necessarily explore the 
backsides of the valley filters. This means that the dislo- 
cation Berry phase gets scrambled, as in the case of ran- 
dom intervalley scattering, and there is no simple rule for 
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FIG. 4: The distributions of the disorder averaged amplitude 
of the harmonics of the magneto-oscillations for the ring de- 
vice of Fig. [T] at zero temperature, with the valley polarizers 
switched off. The black squares indicate the response in the 
absence of a dislocation, and the red stars show what occurs 
in the presence of a d = i dislocation and no intervalley scat- 
tering: the amplitudes of the fundamental and first harmonic 
are halved and their signs reversed (see main text). The tri- 
angles indicate the evolution when the amount of intervalley 
scattering (parametrized by the value of |7p — a? , the prob- 
ability of scattering between valleys on a ring arm traversal, 
expressed in percents) in the arms is increased. 



disentangling the dislocation from the non-topological 
random disorder. At the same time, G(<i>) will be even 
under all circumstances, since there is an infinity of long 
paths in both valleys, and Biittiker's theorem is obviously 
apphcable to this case. 



V. MODELING THE DECOHERENCE AT 
FINITE TEMPERATURE 

We can now turn to the puzzle announced in the In- 
troduction: what happens in our device at finite tem- 
peratures? It appears that our device might represent a 
particular challenge to the incomplete understanding of 
the relation between the coherent quantum transport at 
short scales and classical transport at macroscopic scales 
that is characteristic for any system at a finite temper- 
ature. The sharpest way to express these matters is by 
realizing that at sufficiently large length- and time-scales 
any electron system will be governed by the same hydro- 
dynamical principles as the classical electron plasma of 
the high temperature limit. In contrast to the zero tem- 
perature quantum case, this classical transport is dis- 
sipative and for a Fermi-liquid the dissipation mecha- 
nisms seem well understood; they are the usual electron- 
electron and electron-phonon scattering lore. One can 
just take the Kubo formalism from the textbooks^ and 
compute the diagrams. The problem is that such a com- 
putation becomes unmanageable for a device problem 
such as ours. 

The argument presented in the introduction for the 
unevenness of the magnetoconductance at finite temper- 



ature rests implicitly on the Feynman path^ intuition. 
In the first subsection we will analyze this in more de- 
tail, discovering that the argument actually rests on an 
uncontrolled assumption: to find out what happens with 
the quantum interferences at finite temperature one just 
sums over worldlines up to a maximal length equal to 
the phase coherence length, assuming that the remain- 
der merely contributes to the incoherent current. In this 
way, when worldlines become "too long", they are as- 
sumed to just disappear. In reality these of course do 
not disappear but they turn into the self-energy graphs 
coming from the quasiparticle interactions - the "Kubo 
brick wall" . With this assumption, we obtain the fi- 
nite temperature uneven magnetoconductance (Fig. [5]), 
which becomes even at zero temperature as it should 
(Section |VX|. 

Although it is far from obvious why the cutting of 
world lines approach to dephasing can lead to faulty con- 
clusions regarding the "quantum arrow of time" , prece- 
dents exist where the intuition based on Feynman paths 
turned out to be misleadingii^ It is a standard practice 
in mesoscopic physics to use the scattering matrix ap- 
proach also at finite temperatures, and to account for 
the effects of dephasing using the voltage probe method 
invented by BiittikeriSifSection iy B[). D espite its simplic- 
ity and track record (e.g., Refs. 18 22,23 24,^25), and the 
fact that by construction it respects the basic symmetries 
of quantum scattering, it is surely not a divine solution. 
A problem of principle is of course that this language, 
describing interfering quantum mechanical waves, is not 
quite the preferred way to describe finite temperature 
dissipative fiows of classical hydrodynamics, where the 
quantum unitarity condition is replaced by the weaker 
current conservation demand. In fact, the Biittiker de- 
phasing reservoirs model the effects of inelastic scattering 
by an effective elastic scattering!^ 

For relaxational, classical hydrodynamics, time is at 
the heart of matter. Dealing with a problem like ours, 
where there are subtle complications associated with 
time, can the standard approach be trustworthy in the 
cross-over regime? We favor experimental advances in 
this regime. As we will show in Section IVBl the Biittiker 
construction insists that the magnetoconductance should 
stay even in all circumstances, even when imposing differ- 
ent decoherence rates for the two valleys, as is generally 
expected of this formalism. On the other hand, in the 
high temperature regime the transport turns classical, 
and the expectation of evenness, observed in the large 
body of existing experiments, is theoretically supported 
if given that microscopic reversibility can be viewed as 
certain assumptions on the classical fiuctuation correla- 
tionsjii^i^ 



A. The Feynman path approach 

We describe here the Feynman path approach explic- 
itly."*^ We ignore intervalley scattering of any kind (7 = 
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FIG. 5: The magnetoconductance oscillations G("l>) as a func- 
tion of applied flux $, calculated using the "truncated Feyn- 
man path method" discussed in Section IV Al for the device 
of Fig. [T] with a dislocation of class d — ^. We show the re- 
sults for phase coherence lengths = 3, 7, oo in units of the 
ring arm length, finding that the extremum shifts from ~ i 
the flux quantum value at high temperatures, to the origin at 
zero temperature. The thin dot-dashed line shows the result 
without a dislocation. The inset shows the range of "disor- 
der" (phase (j> of Section IIV|) dependent phase shifts of the 
fundamental period) harmonic of G($) as a function of 
Ly,. Symbols show the average over the "disorder" phase. 

after Eq. ([4])), except at the polarizers, and focus on the 
regime where the phase coherence length is of order of 
the device dimension L. The conductance is proportional 
to the electron transmission probability S, expressed in 
terms of Feynman amplitudes Aa as 

J:^\Ai+A2 + ...+An(lJ^ + B, (6) 

where we assume that in the coherent part only paths 
with a length not exceeding L^p are to be included. 
The longer paths contribute incoherently to the current 
through the term B, i.e. they do not produce interference 
terms responsible for the Aharonov-Bohm oscillations. 
This is the core of the dephasing model of this subsec- 
tion. Let us first discuss the qualitative picture. As al- 
ready explained, the perfect valley polarizer acts by being 
fully transparent to, say, a mode propagating from 
left to right, and a K- mode propagating in the oppo- 
site direction. But microscopic time reversal invariance 
in combination with charge conservation implies that an 
incoming A"_ mode moving to the right is fully reflected 
into a iir_|_ mode moving to the left, and vice versa (in- 
set Fig. [T]). Let us now consider the shortest possible 
paths that can give rise to interference in the presence 
of a dislocation. For a current flowing from left to right, 
the valley polarizers ensure that it is entirely carried by 
modes. The current in the "lower" arm has to tra- 
verse the Volterra cut acquiring the phase jump while in 
the "upper" arm it is unaffected (see Fig. [1]) , with the net 
result that the transmission amplitude picks up the dislo- 
cation pseudoflux of ^, which can in turn be compensated 



by an external field. Repeating the argument for a cur- 
rent from right to left (dashed lines in Fig. [1]), one ends 
up with a shift of — ^ of a flux quantum. The conclusion 
is that the extremum of the magneto-oscillations shifts 
away from its position at zero external flux, thus violat- 
ing Biittiker's theorem. The effect is due to the finite 
temperature and the implicitly dissipative measurement 
setup, and it is present even though we consider the sys- 
tem very close to equilibrium. 

In the present context based on formula ([6|) , the micro- 
scopic time reversal symmetry, which recovers Biittiker's 
law, is associated with the requirement that the perfect 
valley polarizer is a unitary intervalley scatterer for the 
electrons coming in with the wrong polarization. At the 
"high" temperatures discussed in the previous paragraph, 
the phase coherent electron encounters the polarizers at 
most once; it is transmitted with no opportunity to ex- 
plore the "backside" of the polarizers. However, as the 
temperature is decreased one has to take into account 
longer and longer paths. A typical path in this ideal 
device is of the kind that, say, a if+ particle having 
traveled from left to right in the upper and lower arms, 
travels further via the upper arm after getting scattered 
to the valley at the backside of the left polarizer. 
Such long paths destroy the valley quantization and the 
phase associated with the dislocation pseudoflux gets av- 
eraged away. At zero temperature paths of arbitrarily 
long length dominate and the extremum of the magne- 
tooscillation obeys Biittiker's theorem. 

To find out what exactly happens in this model as 
a function of decreasing temperature (and in non-ideal 
devices) , we computed the magnetoconductance by sum- 
ming all Feynman paths with a length bounded by L^p, 
i.e. the part of S without B in Eq. ([6|). The results 
are shown in Fig. [5l as a function of the two parame- 
ters, the dephasing length L^, and the disorder phase (f> 
of Eq. Unsurprisingly, we find a smooth evolution 
where the extremum shifts from a flux « ^ at "high" 
temperature, back to the origin as the phase coherence 
length increases. Only at precisely zero temperature is 
Biittiker's law recovered, since at any finite temperature 
the sum is always "dominated" by the short paths, and 
for this reason the effect seems quite robust. 

The details of the computation go as follows: given the 
finite coherence length measured in ring arm lengths, 
the sum over Feynman paths limited by is performed 
elegantly by a simple trick. We weight the scattering 
matrices of the (single transversal mode) ring arms with 
an auxiliary variable a, essentially scaling t — > at in 
Eq. Then the total coherent transmission ampli- 

tude A(a, $) is calculated exactly in the way described 
in Section llVi by considering the scattering matrices of 
ring arms and polarizers that connect the various in and 
outgoing amplitudes in both valleys, and solving for the 
outgoing amplitudes.— The advantage of doing things 
this way comes from the fact that every Feynman path 
amplitude is exactly a product of scattering matrix ele- 
ments of ring arms, polarizers, and the terminals, which 
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are accumulated as the path is followed from start to 
end4^ Crucially, this implies that the amplitudes of the 
paths having length of n ring arms (traversing an arm 
n times) will pick up the factor a", since a single factor 
a is associated with every pass through an arm. The 
sum of Feynman amplitudes ^'^"^(<i>), corresponding to 
paths of length up to n ring arms, represents the part 
of the total amplitude ^(a, $), where a appears mul- 
tiplied by itself not more than n times. This can be 
obtained by using a truncated Taylor expansion in the 
variable a, because it is an expansion in terms of powers 
of a, exactly what is needed. It follows that = 

, where we used the 
This transmission 



-«M„($) 



nil da'" 



definition Am{^) 

amplitude then gives the conductance associated with 
paths traversing the arms not more than n times through 
the standard relationship G^"),($) ^ 



B. The valley dependent Biittiker dephasing probe 

Let us now turn to the scattering matrix theory 
at finite temperature for our device by employing the 
"Biittiker phantoms" to model the effects of dephasing 
(see Fig. |6]). 

We define T^^' = I^Sp^ 'P, the modulus squared of the 
device scattering matrix elements, where p, q refer to the 
leads (terminals), and a, a' G {+, — } to the propagating 
modes in the leads, such that they represent the proba- 
bility of scattering from mode a' in lead q to mode cr in 
p. It follows that the total current in lead p carried by a 
electrons is 



q,a' 



(7) 



where we use the most general option of having a differ- 
ent chemical potential /i^ for each type cr of electrons, in 
the reservoir connected to the p lead. Such a possibility 
is clearly applicable when we interpret the mathematical 
model as describing a spin system (with two valleys be- 
ing the spin up and down), while for graphene it could 
be less clear what different chemical potentials /i^ 
actually signify. In particular, one could argue that al- 
though the voltage probe is in a sense a mathematical 
construction that enables us to incorporate decoherence 
in the elastic model, it is also an actual component regu- 
larly used in the laboratory, therefore having a physical 
meaning. In the case of graphene, the special voltage 
probe would amount to having different chemical poten- 
tials at the two points in the Brillouin zone, which is 
conceptually conceivable. As will be elaborated below, 
the physical demand for equal dephasing lengths for the 
two valleys leads to ji^ — and removes the problem 
for that situation. In any case, we regard that, conceptu- 
ally, the literal interpretation of the dephasing reservoir 
as a physical entity is not necessary4^ 




FIG. 6: (a) The network of scatterers modeling the 
Aharonov-Bohm device, with dephasing included. Wavy lines 
represent reservoirs and smooth lines represent wires carrying 
the (±) modes. The triangle element and its reservoir belong 
to the Biittiker dephasing probe construction, and are used 
only in Sec. IV Bl the currents Is and ly of Eq. pup are flowing 
in the two leads connecting the triangle to its reservoir. Note 
that we use different chemical potentials (/^j^'s) for the two 
valleys in this reservoir, (b) Labeling of incoming/outgoing 
modes for a generic scatterer, with II representing the column 
etc.; the full/dashed lines depict the +/— modes. 



The Biittiker voltage probe method^^ is based on the 
idea that electrons lose their phase in reservoirs, thus one 
extends the system by introducing N — 2 additional, aux- 
ihary ("phantom") reservoirs (labeled by / e {3,7V}), 
where every one of them is coupled to the device through 
two familiar leads (labeled / and /' at reservoir /), each 
carrying the two ("±") modes, but with the constraint 
that the total current towards a reservoir = 0, i.e. 
the reservoirs will not drain current, but will provide de- 
phasing. The choice of two leads (instead of, e.g. one) is 
just to make possible total decoherence4^ Effectively, one 
solves these N — 2 current constraints (linear equations) 
for the a priori unknown N — 2 auxiliary phantom chem- 
ical potentials ^f, and then eliminates these /i^ in the 
expressions for the currents in the physical leads. Per- 
forming this elimination in the physical current equations 
leads to new, effective transmission coefficients between 
the physical leads which figure in these equations. These 
effective transmission coefficients are then functions of 
the extended system's transmission coefficients between 
the physical leads, as well as the transmission coefficients 
to the phantom leads. To recall the familiar results of 
Ref. [2l|, let us briefly specialize to the simplest, single 
mode, two terminal case, dropping thereby the a index, 
as well as having the simple expression for the conduc- 
tance G = Ti2 implied by Eq. ([7]). Then for example in 
the case of one phantom lead (/ = 3), the above elimi- 
nation procedure yields, 
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where Tg^ is to be understood as the total transmission 
coefficient from lead p to the dephasing reservoir 3, e.g. 
^31 = ^31 + Ts'i, etc. The form of the conductance 
Geff = 712,6// obtained in this way tells us that the 
current divides into a coherent (T'12), and an incoherent 
piece, where the second term is quite suggestive: elec- 
trons starting from lead 2 leave the device to 3 and come 
back to lead 1, while we have to multiply the probabilities 
to obtain the answer — the classical, incoherent way of 
propagation. The amount of decoherence is determined 
by the probability of scattering into the dephasing leads. 
For instance, when electrons leave into the dephaser with 
unit probability, the coherent contribution to the net con- 
ductance vanishes completely, since unitarity of scatter- 
ing requires that T^g — 1, p G {1,2} implies T12 — 
in Eq. ([8]). The bottomline is that the effects of inelas- 
tic scattering are mimicked by a model system of elas- 
tic scatterers with extra leads added, while the current 
constraints become nonlinear in the amplitudes (linear 
in their moduli squared), thereby scrambling the phase 
information. 

An appealing feature of this method is that the ef- 
fective system including the decoherence automatically 
respects the symmetries of the original scattering prob- 
lem;^ in so far as it is encoded in the T matrices. One 
can explicitly check by using formulas of the variety ^ 
(see Ref. 113), that unitarity (sum of elements of any row 
or column of T equals 1), and time reversal symmetry 



-$)), of the starting extended T matrix 



imply precisely the same symmetries for the T^j / matrix. 
These two symmetries are sufficient to derive Biittiker's 
theorem on the evenness of the magnetoconductance for 
a two terminal device 4 

We are now ready to address our graphene device. The 
essential ingredient is the demand that valley currents be 
conserved, such that the Berry phase of the dislocation 
becomes active. The implication is that the phantom 
reservoirs have to respect valley conservation. This is 
at odds with the notion of an equilibrium reservoir that 
would back inject valley currents with equal probability 
regardless the nature of the current it swallows. The 
standard dephasing reservoirs of the Biittiker theory are 
obviously of this equilibrium kind and we have to modify 
the construction to do justice to the conservation law as- 
sociated with the "internal" valley quantum number. We 
first emphasize again that in the standard treatment of 
the single mode case^i with one dephasing probe (labeled 
3), one imposes the hydrodynamical conservation of the 
total current by setting If^^ = 0, which then leads to 
Eq. ([5]) after elimination of the chemical potential /ig. In 
order to allow a maximal current flow in and out of the 
dephasing reservoir, one equips it with two leads labeled 
by / = 3, 3' as we already discussed. Thereby the hydro- 
dynamical current conservation turns into the constraint 



h' - 0, 



(9) 



leads and the two phantom leads can be symbolically 
represented by a triangle as in Fig. [51 

In order to impose the crucial valley current conser- 
vation as well, we now generalize this construction by 
setting 



1^ + 4 



0, 



I-, = 0. 



(10) 



requiring that the dephasing reservoir drains no net cur- 
rent. The scattering matrix connecting the two physical 



In this way we enforce that the decoherence happens in- 
dependently for the two valley currents. 

It is obvious that we have to introduce two chemical 
potentials, and , and use them to enforce the two 
constraints in Eq. pop . Since and /ig" can be used as 
independent parameters, one may conclude that the in- 
trinsic non-equilibration of the conserved valley currents 
is expressed as a non-equilibrium state of the phantom 
reservoir, keeping in mind that in principle this reser- 
voir has no physical existence — it is just a trick to en- 
code that electrons moving through the ring at finite tem- 
perature will dissipate their energy by exciting phonons 
and electron-hole pairs. In summary, the constraints in 
Eq. pn)) are coding for the non-standard ingredient that 
an "internal" (valley) quantum number is conserved. Fi- 
nally, we emphasize that there is an additional freedom in 
the choice of the scattering matrix (5^;) associated with 
the way the dephaser is connected to the ring: the scat- 
tering indicated by the triangle element of Fig. [51 This 
contains the transmission coefficients into the dephaser, 
and thereby controls the degree of decoherence caused 
by the dephaser. The does not mix the valleys, so 
it determines the intravalley dephasing time. Obviously 
in the physical system the decoherence in the two val- 
leys should be the same, leading to constraints on the 
matrix elements discussed in detail below. Given these 
ingredients, the calculations are straightforward and are 
summarized in the Appendix. 

The outcome for magnetoconductance computed from 
the results in the Appendix is as follows. According to 
expectations, we analytically prove that the magnetocon- 
ductance G{^) is even in the flux $, assuming the sym- 
metries (time reversal and unitarity) of the T matrix are 
present. The evenness is thus independent of the values 
of all physical parameters, and persists even when differ- 
ent dephasing lengths are assigned to the two valleys by 
tuning the matrix. Such a situation corresponds to 
a non-equilibrated reservoir, with /i^ 7^ Ma"- Invariably, 
these chemical potentials scale with the physical voltage 
Ml ~ (no ± dependence in physical reservoirs) , consis- 
tent with the linear response regime. 

Let us now analyze the oscillations themselves. At 
zero temperature, when the scattering into the phan- 
tom reservoirs vanishes, the model reduces to the mat- 
ters discussed in Section [TVl The corresponding results 
for the disorder phase averaged amplitude of the funda- 
mental, ^ harmonic, seen as the flrst entry of Fig. [31 are 
shown as the infinite dephasing length (L^ = 00) entry 
in Fig. [71 (black square — in absence of dislocation, red 
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L^[ring arm] 



FIG. 7: The temperature, parametrized through the dephas- 
ing length L^, dependence of the disorder phase averaged 
harmonic amphtude, within the Biittiker dephasing modeL 
The symbols are taken from Fig. |4] — black square denoting 
the ideal ring, red star the dislocated ring with no interval- 
ley scattering, and triangles the dislocated ring with varying 
intervalley scattering lengths. The first, = oo entry repro- 
duces the zero temperature result of Fig. [l] At each separate 
finite dephasing length value, amplitudes are normalized by 
the ideal ring amplitude (black square). The figure then shows 
how the ratio of the dislocated ring harmonic (red star) to the 
ideal ring harmonic varies negligibly from the zero tempera- 
ture value of —1/2. The green dots represent the values of 
the ideal ring un-normalized amplitudes at each dephasing 
length, and show how the oscillations disappear with rising 
temperature. 



star — in presence of dislocation, triangles — with dis- 
location and varying intervalley scattering length). The 
effect of finite temperature is modeled by switching on 
the scattering into the dephasing reservoir, and amounts 
to a gradual decrease of the magnetoconductance oscilla- 
tions that eventually vanish when the dephasing length 
becomes small compared to the device dimensions: the 
green dashed line and green circles of Fig. [7] show the 
overall oscillation amplitude dependence on the dephas- 
ing length L^. The next issue is how the ratio between 
the disorder phase averaged Fourier amplitudes of the 
dislocated and ideal ring evolve with temperature. The 
Fig. [7] shows that this ratio is virtually independent of 
the temperature, and retains the value of — ^ identified 
at zero temperature (Section lIVp . 



VI. CONCLUSIONS 



a more complete theoretical treatment that has the ca- 
pacity to settle these matters. We hope that the consid- 
erations in this paper will motivate the experimentalists 
to focus in on the physics of dislocations in graphene. 
It seems there are no fundamental obstructions to the 
realization of our proposed device, with the possible ex- 
ception that it might appear challenging to keep valley 
currents conserved on reasonable length scales. On the 
other hand, such an experiment still represents a consid- 
erable technical challenge, but the reward is potentiality 
of probing the reach of validity of a familiar law in a novel 
context. 

We do invoke specialties of graphene, but the theme 
is more general. Are there other conserved internal 
quantum numbers that can be utilized for similar pur- 
poses? The transport of spin comes immediately to 
mind, with spin polarization taking the role of valley 
polarization, spin currents^iiS as valley currents, and 
spin-orbit coupling being like intervalley scattering. One 
needs more equipment. It appears that in principle the 
Aharonov-Casher Berry phased associated with an elec- 
trical monopole in the middle of the ring has the potential 
to take the role of the dislocation (see, e.g., Ref. [13), but 
a literal analogue of valley polarizers is less obvious. 

This brings us back to an important byproduct of this 
pursuit: the graphene dislocation with its Berry phase 
that communicates with valley currents in a unique way. 
More speculatively, if "valleytronics" ever gets off the 
ground, and the Casimir-Onsager relations are shown 
to fail in the intermediate regime (however unlikely the 
prospect), the dislocations would have their use as unique 
valleytronic circuit elements measuring in a topologically 
robust ways the valley polarized currents. The equipment 
based on valley filters, that were at the focus of this pa- 
per, might not be the best way to go, and the same objec- 
tion holds for other possible microscopic mechanisms of 
producing valley polarized states An analogy with the 
quantum spin Hall effect^i^i^i^ suggests another alley 
to explore. Topological band insulation rooted in spin- 
orbit coupling goes hand in hand with chiral spin currents 
at the surface, and it is imaginable^° that these can be 
exploited to construct a spin battery. It was recently ar- 
gued that similar topologically protected currents exist at 
the interface between graphene bilayers, where the gap 
associated with AB sublattice breaking changes sigui^ 
These chiral interface states are associated with valley 
currents and one can contemplate to exploit these for 
the purpose of constructing a valley battery. 



Can Casimir-Onsager relations be invalidated because 
leftovers of quantum coherent currents at intermediate 
temperatures cannot equilibrate in a true sense due to 
a conservation law applying to an internal symmetry? 
The special features of the device introduced in this pa- 
per make this provocative question germain. We do not 
claim to have a definitive answer. Within the realm of fi- 
nite temperature quantum transport the issue appears to 
be unresolved and we challenge the readership to devise 
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VII. APPENDIX 

We discuss here the details of the scattering matrix 
calculations which include the valley preserving dephas- 
ing reservoirs, as discussed in Section IV Bl We follow 
in detail the method and interpretation introduced in 
Refs. Pl50|. We eliminate the dephaser, keeping every- 
thing else in the system arbitrary. At this point we 
need only the information that there are two time re- 
versed modes in the system, while the dephaser takes 
the form (fTO|) . and the full T matrix respects the basic 
symmetry requirements, i.e. unitarity and time rever- 
sal. We label the eliminated dephaser by / = 3, keeping 
the label p for all the other leads of the device, some 
of which might refer to other dephasers. By applying 
the expression for the current in an arbitrary lead ([7]) 
to the dephaser (i.e. putting p = 3,3'), we determine 
the phantom potentials /i^ , , which are present in the 
constraint equations ()10p . in terms of the other poten- 
tials. Once these phantom potentials are determined, 
we eliminate them in the expressions for the currents in 
all the other terminals of the device (i.e. p ^ 3,3' in 
Eq. ([7])). These expressions then describe the currents in 
all the remaining terminals, as a function of their associ- 
ated chemical potentials in the form of an effective matrix 
T^f f with elements T^ff pg ■ The general expressions for 
the elements of Te//, and the chemical potentials of the 
eliminated dephasing reservoir are 



translated into the equality J2qf3 a ^'^g'^ ~ ^' ^^"^ ^^^^ 
equality is proved correct directly by using the unitarity 
of the starting T matrix, i.e. ^, T^^ = 1. Time rever- 
sal symmetry of the original matrix is expressed through 
T^^'{<^) = r*p'*(-$), as is impHed by the time rever- 
sal property of the matrix S_, Eq. It is straightfor- 
ward to show that this property also holds for T^ff of 
Eq. pT|) . by checking it for explicit values of a and a', 
with use of the property for T; a simplification comes 
from noting that A(<1>) = A(— $). This finishes our anal- 
ysis of the general case. For the particular case of our 
two terminal device model with one dephasing element 
in the ring arm. Fig. [6l we also evaluated the conduc- 
tance Gphys = J2a(T' numerically, and the results 
are discussed in Section IV Bl We crudely estimate the 
dephasing length L^, used qualitatively in Fig. [3 by the 
simple formula exp{—l/L,^) = 1 — where e is the 
element of the 5^ matrix which describes the scattering 
from the ring arm into the dephasing reservoir. The value 
of 1 — jep is then the probability for the electron not to 
dephase while traversing the arm. By considering each 
ring arm traversal (path of length 1) as an independent 
statistical measurement of this probability, the charac- 
teristic dephasing length follows. This is analogous 
to writing the formula for the half life of an unstable 
particle having the probability jep to decay. 



rrifTfT rp(T(T 

eff,pq ~ ^l^q 



MP" 



(11) 



2 — Tff^Tf"^ + TP-T-'^ 

33 ^ 3q 33 3q 



A^(2-T±+)(2-T33^ 



rp-\ rp 1- 

33 33 ' 



where a denotes the valley opposite to a {a — ^ for a = 
±), and we used the obvious abbreviations for summing 
over the leads (3, 3'e 3), e.g. Tj' = T^/ ^Tg, Tg = 

rpaa' I rpaa' _i_ rpcrn' _i_ rpatr' „j.„ 

One can in principle now proceed to eliminate the next 
dephasing circuit element, using the same procedure out- 
lined above, of course now applied to the effective matrix 
Teff defined in Eq. (|lip . which describes the system at 
this stage. If this is repeated iteratively for all phan- 
toms, the final effective Tphys, which describes scatter- 
ing between the physical leads, is obtained. The general 
implication of this dephasing model, proven in the fol- 
lowing, is that the two basic symmetries (unitarity and 
time reversal) hold for Teff in Eq. ()lip . which in turn 
implies that they will hold also for each T matrix ob- 
tained by successive eliminations of dephasing reservoirs, 
including the sought Tphys at the final stage. If Tphys 
describes a two terminal device, we conclude that the 
assumptions of Biittiker's theorem hold, and the magne- 
toconductance is even. The unitarity of scattering de- 
scribed by Teff, expressed as J2q a' "^eff pq — 1' '"^'^ ^le 



To clarify the structure of the quantities calculated 
above, Eq. ([TT|) . and compare them to the one mode case 
of Ref. 4, we now specialize to the relevant two terminal 
(labeled interchangeably by L, R for "left" and "right" , 
and 1, 2, respectively) device with one additional de- 
phasing reservoir 3. It will then be useful to introduce 
quantities which group T matrix elements according to 
their meaning, so we define: the total probability of scat- 
tering the electron from the physical lead p into a a 
electron exiting into the dephasing reservoir by Sf^ = 
Tip^ + ^ap" -I- Ty+ + TgTp", the probabiHty of scattering 
an incoming a electron from the dephasing reservoir into 



the p physical lead by = ^ij'^ 



rp-\-(7 I rp — a 



as well as the obvious quantities S^, — S^t + S^, , 



aa clcr 



Slj^, and = 



ij c 



ij c 



Sqo- A physi- 



cal quantity of relevance is the amount of energy dissi- 
pated by the dephaser, which acts as an inelastic scat- 
terer,— by the cr current. This describes the amount of 
dephasing of the a current, and we express it through 
a dimensionless function rf as = j^ri'^{^L ~ Mfl,)^- 
It reflects the fact that the electrons exchanged between 
the physical and dephasing reservoir are injected at dif- 
ferent chemical potentials. If one also defines dimen- 
sionless functions x"^ to describe the phantom chemical 
potentials by fi^ = fin + x^it^-L — fJ-n), and introduce 
one more auxiliary combination of matrix elements by 
Cp„ = S^Sf + TgSf, then for the two terminal device 
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it can be written that 



X 



2+ 



r^- = S}^ + Tg{x^f - (2 - ^g^)(x")^ 

Teff,21 = T2I + Sl^X^ + Sl^X^ = Teff,12- 



(12) 



Physically, one expects the degree of decoherence to be 
the same for the two valleys, but we emphasize again that 
this is not needed for the evenness of magnetoconduc- 
tance, since no such assumption was made in its proof 
at the beginning of the Appendix. This demand reads 
= W~ , and in general puts constraints on the ma- 
trix elements. Although it is not trivial to extract the 
simplest condition equivalent to this demand, we note 
that in the case when 5*^^+ — S}^ = ^S} = S}, and 
5f+ = Sf- = ^Sf = Si (the use of underline S symbols 
here should of course not confuse with the same symbol 
for scattering matrices), which corresponds to saying that 
the probability of the electron coming from the first lead 
to scatter into the dephasing reservoir as + or — electron 
is the same, and this statement holds also for the second 



lead, significant simplifications occur. Namely, if we put 
K/o = ^i/o = l^i/o = ^i/o^ the quantities in Eq. 
become. 



IT =7] = 



e//,21 



21 



C2 cl 



(13) 



The two valleys contribute equally to the incoherent part 
of the transmission in the last equality. Finally, we note 
that if additionally to the above assumption, which leads 
to = W~ , we assume a vanishing magnetic field, 
the time reversal condition 5f ($) — S'o(— <&) becomes 
Sf — Sg, and we see from (IT^ that the dimcnsionless 
dissipated energy becomes equal to the incoherent trans- 
mission, i.e. 2ri = T^ff ^21 — T21. The factor 2 accounts 
for two modes; this recovers the result of Ref. [2l|, which 
considered a single mode and zero magnetic field situa- 
tion. 
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